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(57)Abstract: 

PROBLEM TO BE SOLVED: To improve the accuracy in operating 
the amount of the air charged in a cylinder in transition. 
SOLUTION: An opening command value is determined 
corresponding to an accelarator operation amount or the like during 
the operation of an engine, and the output timing of the opening 
command value is delayed by a predetermined delay time. An 
estimated variation of a throttle opening is operated by an 
electronic throttle valve on the basis of the opening command 
value before delay, the estimated variation is added to the present 
throttle opening (the output of the throttle opening sensor) to 
determine an estimated throttle opening at an intake valve closing 
timing. Then the estimated amount of the air charged in the 
cylinder is temporarily operated by the intake system model on the 
basis of the estimated throttle opening, and the obtained amount is 
differentiated and integrated to operate the estimated variation of 
the amount of the air charged in the cylinder. The estimated 
variation is added to the amount of the air charged in a base 
cylinder operated by a base intake system model to determine the 
final estimated amount of the air charged in the cylinder (the 
amount of the air in the cylinder determined at the intake valve closing timing). 
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* NOTICES * 

iTPO and NCI PI are not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 
2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



CLAIMS 



[Claim(s)] 

[Claim 1] In the internal combustion engine having the electronic throttle system which drives a throttle valve 
with a throttle actuator and controls a throttle opening An opening command value operation means to calculate 
an opening command value based on an accelerator control input etc., A delay means to delay the timing which 
outputs the opening command value calculated with said opening command value operation means to said 
throttle actuator, A throttle opening prediction means to predict a subsequent throttle opening before the delay 
output of this opening command value based on the opening command value before making it delayed with said 
delay means, and said electronic throttle system response delay property, A restoration air content prediction 
means in a cylinder to predict the restoration air content in a cylinder based on the throttle opening predicted 
with said throttle opening prediction means, The control unit of the internal combustion engine characterized by 
having a fuel-oil-consumption operation means to calculate fuel oil consumption based on the restoration air 
content in a cylinder predicted with said restoration air content prediction means in a cylinder. 
[Claim 2] Said restoration air content prediction means in a cylinder is the control unit of the internal 
combustion engine according to claim 1 characterized by predicting the variation of the restoration air content 
in a cylinder to intake valve close timing based on the throttle opening predicted with said throttle opening 
prediction means, adding to the restoration air content in a base cylinder which calculated this variation based 
on the current operation parameter, and predicting the restoration air content in a cylinder. 
[Claim 3] Said restoration air content prediction means in a cylinder is the control unit of the internal 
combustion engine according to claim 2 characterized by to predict the variation of the restoration air content in 
a cylinder to intake valve close timing using the inhalation-of-air system model which applied the law of mass 
conservation to the inhalation air which considers that throttle opening which inhalation air passes is an orifice, 
and flows a throttle passage air content and a throttle down-stream path in integrating the variation of the output 
of this inhalation-of-air system model to intake valve close timing. 

[Claim 4] The formula which calculates a throttle passage air content with said inhalation-of-air system model 
is [Equation 1]. 
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In case it is alike, and is set up and said restoration air content prediction means in a cylinder calculates said 
throttle passage air content, f (Pm/Pa) is Pm/Pa. It is the control unit of the internal combustion engine 
according to claim 3 characterized by computing from the table made into a parameter and computing mu-A 
from the table which makes a throttle opening a parameter. 

[Claim 5] The table of f (Pm/Pa) used for the operation of said throttle passage air content It is set as an f 
(Pm/Pa) = negative value at the time of f(Pm/Pa) =0Pm/Pa >1 at the time of f(Pm/Pa) = forward value Pm/Pa 
=1 at the time of Pm/Pa <1 . Said restoration air content prediction means in a cylinder is the control unit of the 
internal combustion engine according to claim 4 characterized by having a means to equalize the operation 
value of said inhalation-of-air system model. 

[Claim 6] Said delay means is the control unit of the internal combustion engine according to claim 1 to 5 
characterized by setting the time delay of an opening command value as the time amount which deducted the 
dead time of said electronic throttle system from the time amount from the operation timing of the fuel oil 
consumption of a certain cylinder to the intake valve close timing of the cylinder concerned. 
[Claim 7] Said delay means is the control unit of the internal combustion engine according to claim 1 to 6 
characterized by outputting without delaying said opening command value when the time amount from the 
operation timing of the fuel oil consumption of a certain cylinder to the intake valve close timing of the cylinder 
concerned becomes shorter than the dead time of said electronic throttle system. 

[Claim 8] Said delay means is the control unit of the internal combustion engine according to claim 1 to 7 
characterized by outputting without delaying said opening command value when an automatic transmission 
corresponds to either of the neutral conditions at the inside of the predetermined time immediately after start up, 
and the time of idle operation at the time of start up. 

[Claim 9] Said throttle opening prediction means is the control unit of the internal combustion engine according 
to claim 1 to 8 characterized by predicting the throttle opening after an opening command value delay output 
using the electronic throttle model containing the lag unit and rate limiter more than primary [ which considers 
the opening command value before making it delayed with said delay means as an input ]. 
[Claim 10] Said throttle opening prediction means is the control unit of the internal combustion engine 
according to claim 1 to 9 characterized by predicting the variation of the throttle opening to intake valve close 
timing using said electronic throttle model, adding this variation to a current throttle opening, and predicting the 
throttle opening of intake valve close timing. 

[Claim 11] Said fuel-oil-consumption operation means is the control unit of the internal combustion engine 
according to claim 1 to 10 characterized by having a means to amend fuel oil consumption according to 
operational status, and switching the correction factor to fuel oil consumption by the time of the load effect by 
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accelerator actuation, and the case of being other. 

[Claim 12] In the internal combustion engine having the electronic throttle system which drives a throttle valve 
with a throttle actuator and controls a throttle opening An opening command value operation means to calculate 
an opening command value based on an accelerator control input etc., A throttle opening prediction means to 
predict the throttle opening of intake valve close timing based on the opening command value calculated with 
said opening command value operation means, and said electronic throttle system response delay property, A 
restoration air content prediction means in a cylinder to predict the restoration air content in a cylinder based on 
the throttle opening predicted with said throttle opening prediction means, The control unit of the internal 
combustion engine characterized by having a fuel-oil-consumption operation means to calculate fuel oil 
consumption based on the restoration air content in a cylinder predicted with said restoration air content 
prediction means in a cylinder. 

[Claim 13] Said restoration air content prediction means in a cylinder is the control unit of the internal 
combustion engine according to claim 12 characterized by predicting the variation of the restoration air content 
in a cylinder to intake valve close timing based on the throttle opening predicted with said throttle opening 
prediction means, adding to the restoration air content in a base cylinder which calculated this variation based 
on the current operation parameter, and predicting the restoration air content in a cylinder. 

[Claim 14] A restoration air content operation means in a base cylinder to calculate the restoration air content in 
a base cylinder based on a current operation parameter, The inhalation-of-air system model which applied the 
law of mass conservation to the inhalation air which considers that throttle opening which inhalation air passes 
is an orifice, and flows a throttle passage air content and a throttle down-stream path is used. A variation 
prediction means to predict the variation of the restoration air content in a cylinder to intake valve close timing 
based on the variation of the output of this inhalation-of-air system model, A restoration air content prediction 
means in a cylinder to add the variation predicted with said variation prediction means to the restoration air 
content in a base cylinder calculated with said restoration air content operation means in a base cylinder, and to 
predict the restoration air content in a cylinder, The control unit of the internal combustion engine characterized 
by having a fuel-oil-consumption operation means to calculate fuel oil consumption based on the restoration air 
content in a cylinder predicted with said restoration air content prediction means in a cylinder. 
[Claim 1 5] An intake-air-flow detection means to detect the flow rate of the inhalation air which flows an 
internal combustion engine's inhalation-of-air path, The inhalation-of-air system model which simulated the 
behavior of inhalation air until the inhalation air which passed the throttle valve flows in a cylinder is used. It 
has an operation means to input the output of said intake-air- flow detection means into this inhalation-of-air 
system model, and to calculate the restoration air content in a cylinder which is the output of this inhalation-of- 
air system model. The time constant of said inhalation-of-air system model The control unit of the internal 
combustion engine characterized by setting air content change as a small value so that a actual twist may also 
appear at an early stage. 

[Claim 16] A means to presume the current restoration air content in a cylinder based on a current throttle 
opening, A throttle opening prediction means to predict a future throttle opening, and a means to predict the 
future restoration air content in a cylinder based on said future throttle opening, A means to add to the 
restoration air content in a base cylinder which calculated the deflection of said future restoration air content in 
a cylinder, and said current restoration air content in a cylinder based on the current operation parameter, and to 
calculate the final restoration air content in a prediction cylinder, The control unit of the internal combustion 
engine characterized by having a fuel-oil-consumption operation means to calculate fuel oil consumption based 
on said final restoration air content in a prediction cylinder. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] This invention relates to the control unit of the internal combustion engine which 
improved the operation approach of an internal combustion engine's restoration air content in a cylinder. 
[0002] 

[Description of the Prior Art] In order to correspond to the emission gas clarification regulation which becomes 
severer future still, it is necessary to high-degree-of-accuracy-ize Air Fuel Ratio Control (fuel-injection control), 
to calculate the air content (restoration air content in a cylinder) with which it fills up in an engine cylinder with 
a sufficient precision for that purpose, and to set up the proper fuel oil consumption corresponding to the 
restoration air content in a cylinder. The operation approach of current and the restoration air content in a 
cylinder currently generally used detects an intake air flow with the air flow meter installed in the upstream of a 
throttle valve, and is divided roughly into the method (mass flow method) which calculates the restoration air 
content in a cylinder from the detection value, and the method (speed density method) which detects an intake 
pressure by the intake-pressure sensor, and calculates the restoration air content in a cylinder from an intake 
pressure and an engine speed. 
[0003] 

[Problem(s) to be Solved by the Invention] By the way, the timing which the restoration air content in a 
cylinder decides is intake valve close timing which an intake stroke ends, and the timing which computes fuel 
oil consumption is a front [ timing / intake valve close ] (it is because it is necessary to perform fuel injection 
before intake valve close timing to make an injection fuel inhale in a cylinder). Even if it computes the 
restoration air content in a cylinder by which approach of the conventional mass flow method mentioned above 
and a speed density method, therefore, to a transient The restoration air content in a cylinder changes from the 
calculation timing of fuel oil consumption before intake valve close timing (definite timing of the restoration air 
content in a cylinder). Consequently, the ratio (air-fuel ratio) of the actual restoration air content in a cylinder 
and the fuel quantity which flows in a cylinder shifts from a target air- fuel ratio, and there is a fault that the Air 
Fuel Ratio Control precision of a transient worsens. 

[0004] This invention is made in consideration of such a situation, therefore the object is in offering the control 
unit of the internal combustion engine which can raise the Air Fuel Ratio Control precision of a transient. 
[0005] 

[Means for Solving the Problem] In order to attain the above-mentioned object, this invention predicts the 
throttle opening of intake valve close timing (definite timing of the restoration air content in a cylinder), 
predicts the restoration air content in a cylinder based on the prediction throttle opening, and calculates fuel oil 
consumption based on the restoration air content in a prediction cylinder. Here, the change of the restoration air 
content in a cylinder of the reason using a throttle opening as a parameter which predicts the restoration air 
content in a cylinder is mainly because change of a throttle opening is generated as the origin and change of the 
restoration air content in a cylinder can be predicted with a sufficient response from change of a throttle opening 
to a transient. 

[0006] It becomes possible to predict the throttle opening of intake valve close timing (definite timing of the 
restoration air content in a cylinder) by delaying moderately the timing which outputs an opening command 
value to a throttle actuator with a delay means in the internal combustion engine having the electronic throttle 
system which drives a throttle valve with a throttle actuator and controls a throttle opening like claim 1 . Under 

http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje 1 0/1 7/2005 



JP 3 2002-201998,A [DETAILED DESCRIPTION] 



Page 2 of 17 



the present circumstances, what is necessary is just to predict a subsequent throttle opening before the delay 
output of this opening command value based on the opening command value before making it delayed with a 
delay means, and an electronic throttle system response delay property, since response delay (dead time) exists 
in actuation of an electronic throttle system. It can become possible to predict the throttle opening of intake 
valve close timing with a sufficient precision, and to predict the restoration air content in a cylinder with a 
sufficient precision from the prediction throttle opening by this, and the Air Fuel Ratio Control precision of a 
transient can be raised. 

[0007] By the way, the approach of predicting the restoration air content in a cylinder from a prediction throttle 
opening has the advantage that the responsibility of a transient is good, but on the other hand it has the 
inclination for the forecast at the time of a stationary to shift from a actual value according to dispersion in an 
electronic throttle system, aging, a service condition, etc. Moreover, at the time of a stationary, since the 
restoration air content in a cylinder does not change, the restoration air content in a cylinder calculated based on 
current operation parameters (an intake air flow, intake pressure, etc.) is in agreement with the restoration air 
content in a cylinder decided to subsequent intake valve close timing. 

[0008] Then, it is good to predict the variation (a part for transient change) of the restoration air content in a 
cylinder to intake valve close timing based on a prediction throttle opening, to add to the restoration air content 
in a base cylinder which calculated this variation based on the current operation parameter, and to predict the 
restoration air content in a cylinder like claim 2. If it does in this way, it will become possible to predict the 
restoration air content in a cylinder with a sufficient precision by both the time of a stationary, and the transient. 

[0009] Furthermore, it is good to predict the variation of the restoration air content in a cylinder to intake valve 
close timing using the inhalation-of-air system model which applied the law of mass conservation to the 
inhalation air which considers like claim 3 that throttle opening which inhalation air passes is an orifice, and 
flows a throttle passage air content and a throttle down-stream path in integrating the variation of the output of 
this inhalation-of-air system model to intake valve close timing. By using such an inhalation-of-air system 
model, the variation of the restoration air content in a cylinder to intake valve close timing can be predicted with 
a sufficient precision by comparatively easy data processing. 

[0010] In this case, the formula which calculates a throttle passage air content with an inhalation-of-air system 
model is good like claim 4 to use a degree type. 
[0011] 
[Equation 2] 

Gin = ti - A • Pa • f (Pm/Pa) 
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[0012] In case a throttle passage air content is calculated, f (Pm/Pa) is Pm/Pa. It computes from the table made 
http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web cgi ejje 1 0/1 7/2005 



JP,2002-201998,A [DETAILED DESCRIPTION] 



Page 3 of 17 



into a parameter, and you may make it compute mu-A from the table which makes a throttle opening a 
parameter. Thereby, data processing of an inhalation-of-air system model type becomes very easy. 
[0013] Like claim 5, furthermore, the table of f (Pm/Pa) It is good to set it as an f(Pm/Pa) = negative value at 
the time of f(Pm/Pa) =0Pm/Pa >1 at the time of f(Pm/Pa) = forward value Pm/Pa =1 at the time of Pm/Pa <1, 
and to equalize the operation value of an inhalation-of-air system model. If f (Pm/Pa) is set to f(Pm/Pa) = 0 at 
the time of Pm/Pa >1 although it does not become a negative value physically so that it may mention later, they 
are Pm/Pa. There is an inclination which the operation value of an inhalation-of-air system model vibrates, and 
hunting generates at the time of heavy load operation changed in the one neighborhood. This cause is Pm/Pa. 
They are [ that the rate of change of f (Pm/Pa) becomes large in the field of the one neighborhood, and ] 
operation top Pm / Pa at the time of heavy load operation. Since f (Pm/Pa) is guarded by 0 whenever it becomes 
one or more, it is because change of f at the time of heavy load operation (Pm/Pa) becomes irregular. 
[0014] As this cure, it is f(Pm/Pa) = negative a value, then Pm/Pa like claim 5 at the time of Pm/Pa >1. At the 
time of heavy load operation changed in the one neighborhood, change of f (Pm/Pa) will become regular. For 
this reason, by equalizing the operation value of an inhalation-of-air system model, the operation value of the 
inhalation-of-air system model at the time of heavy load operation can be stabilized, and hunting can be 
prevented. 

[0015] moreover, time delay Tdly of the opening command value delayed with a delay means like claim 6 Time 
amount Tinj from the operation timing (prediction timing of the restoration air content in a cylinder) of the fuel 
oil consumption of a certain cylinder to the intake valve close timing of the cylinder concerned from — it is 
good to set it as the time amount (Tdly =Tinj-Tth) which deducted the dead time Tth of an electronic throttle 
system. It is the time delay Tdly of an opening command value so that the throttle opening of intake valve close 
timing may be made in agreement with a prediction throttle opening, if it does in this way. It can set up and the 
operation of a prediction throttle opening becomes easy. 

[0016] In this case, the dead time Tth of an electronic throttle system is the time amount Tinj from the operation 
timing (prediction timing of the restoration air content in a cylinder) of fuel oil consumption to intake valve 
close timing, although it does not change even if a throttle actuation rate changes. It becomes so short that an 
engine speed becomes a high revolution. For this reason, at the time of a high revolution, it is the time amount 
Tinj from the operation timing of fuel oil consumption to intake valve close timing. There is a possibility of 
becoming shorter than the dead time Tth of an electronic throttle system. 

[0017] This point is taken into consideration and it is the time amount Tinj from the operation timing of fuel oil 
consumption to intake valve close timing like claim 7. When becoming shorter than the dead time Tth of an 
electronic throttle system, it is good to make it output without delaying an opening command value. If it does in 
this way, while not performing useless throttle delay control at the time of a high revolution, the throttle 
responsibility at the time of a high revolution can also be raised. 

[0018] Moreover, like claim 8, when an automatic transmission corresponds to either of the neutral conditions 
in the predetermined time immediately after start up at the time of idle operation at the time of start up, it is 
good to make it output without delaying an opening command value. The time of start up, and immediately 
after start up, essentially, since the engine revolution is unstable, when throttle delay control by which an 
opening command value is delayed is performed, there is a possibility that engine revolution fluctuation may 
become still larger. Moreover, at the time of idle operation, in order for idle rotational-speed control (ISC) to 
operate and to carry out feedback control of the idle rotational speed, when throttle delay control is performed, 
there is a possibility that idle rotational-speed control may interfere with throttle delay control, and an idle 
revolution may become instability. Moreover, since an operator may perform racing (an engine carrying out a 
fore blow) when an automatic transmission is in a neutral condition, when it is in a neutral condition and 
throttle delay control is performed, the start of an engine speed is overdue at the time of racing activation, and 
there is a possibility that accelerator responsibility and acceleration nature may impress an operator with it 
being bad. 

[0019] Therefore, if it is made not to perform throttle delay control like claim 8 when it is the operational status 
in which the adverse effect by throttle delay control appears (at the time [ At the time / At the time of start up / 
of idle operation ] of a neutral), the adverse effect by throttle delay control can be lost. 

[0020] Moreover, in case the throttle opening after the delay output of this opening command value is predicted 
like claim 9 using the opening command value before making it delayed with a delay means, it is good to 
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predict a throttle opening using the electronic throttle model containing the lag unit and rate limiter more than 
primary [ which considers the opening command value before making it delayed as an input ]. Although it is 
generally difficult to carry out physical modeling of this at a precision since an electronic throttle system has 
complicated structure By simulating an electronic throttle system response delay property with the lag unit more 
than primary, and simulating the marginal property of the actuation rate of a throttle valve by the rate limiter 
Even if it can build the easy electronic throttle model of data processing and does not high-performance-ize 
especially CPU of an electronic throttle system, the prediction operation of a throttle opening is attained. 
[0021] Moreover, a throttle opening forecast may shift from a actual value according to dispersion in an 
electronic throttle system, aging, a service condition, etc. Then, the variation of the throttle opening to intake 
valve close timing is predicted using an electronic throttle model, and this variation is added to a current throttle 
opening, and you may make it predict the throttle opening of intake valve close timing like claim 10. If it does 
in this way, it will become possible to lessen the throttle opening prediction error by the above-mentioned 
cause, and to predict a throttle opening with a sufficient precision. 

[0022] Moreover, in case fuel oil consumption is amended according to operational status, you may make it 
switch the correction factor to fuel oil consumption like claim 1 1 by the time of the load effect by accelerator 
actuation, and the case of being other. That is, since this invention can predict the restoration air content in a 
cylinder with a sufficient precision to the load effect by accelerator actuation, it can lessen amendment to fuel 
oil consumption. However, since the load effect by the time of shifting an automatic transmission from a neutral 
range to a drive range, power steering, a brake, an air-conditioner, etc. cannot be predicted from an accelerator 
control input, it is desirable to make [ many ] amendment to fuel oil consumption to these load effects. 
[0023] Therefore, if the correction factor to fuel oil consumption is switched like claim 1 1 by the time of the 
load effect by accelerator actuation, and the case of being other, amendment of fuel oil consumption can be 
rationalized according to the cause of a load effect. 

[0024] Although invention concerning claims 1-11 explained above was made to perform throttle delay control, 
based on the opening command value and the electronic throttle system response delay property which were 
calculated based on the accelerator control input etc., the throttle opening of intake valve close timing is 
predicted, and you may make it predict the restoration air content in a cylinder like claim 12 based on this 
throttle opening, without performing throttle delay control. Even if such, it can become possible to predict a 
throttle opening using the dead time of an electronic throttle system, and to predict the restoration air content in 
a cylinder with a sufficient precision from the prediction throttle opening, and the Air Fuel Ratio Control 
precision of a transient can be raised. 

[0025] It is good to add to the restoration air content in a base cylinder which predicted the variation (a part for 
transient change) of the restoration air content in a cylinder to intake valve close timing like claim 1 3 based on 
the prediction throttle opening, and calculated this variation also in this case based on the current operation 
parameter, and to predict the restoration air content in a cylinder. If it does in this way, it will become possible 
to predict the restoration air content in a cylinder with a sufficient precision by both the time of a stationary, and 
the transient. 

[0026] Although claims 1-13 explained above apply this invention to an internal combustion engine with an 
electronic throttle system In the case of the mechanical-cable-type throttle system which interlocks a throttle 
opening with accelerator actuation mechanically Like claim 14, while calculating the restoration air content in a 
base cylinder based on a current operation parameter The inhalation-of-air system model which considers that 
throttle opening is an orifice and calculates a throttle passage air content from a current throttle opening etc. is 
used. Based on the variation of the output of this inhalation-of-air system model, the variation of the restoration 
air content in a cylinder to intake valve close timing is predicted, this variation is added to the restoration air 
content in a base cylinder, the restoration air content in a cylinder is predicted, and you may make it calculate 
fuel oil consumption based on this restoration air content in a cylinder. Also in the case of a mechanical-cable- 
type throttle system, if it does in this way, the operation precision of the restoration air content in a cylinder can 
be raised conventionally, and the Air Fuel Ratio Control precision of a transient can be raised. 
[0027] Moreover, although the inhalation-of-air system model which calculates the restoration air content in a 
cylinder from a throttle opening was used in above-mentioned claims 1-14, when using the inhalation-of-air 
system model which calculates the restoration air content in a cylinder from the output (inhalation air content) 
of an intake-air- flow detection means (air flow meter), it is good like claim 1 5 for air content change to set the 
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time constant of an inhalation-of-air system model as a small value also for a actual twist so that it may appear 
at an early stage. Thus, if the time constant of an inhalation-of-air system model is made small, the effectiveness 
as predicting the future restoration air content in a cylinder that a actual twist also has it since change of the 
restoration air content in a cylinder calculated with the inhalation-of-air system model appears at an early stage 
will be acquired. [ same ] Thereby, the operation precision of the restoration air content in a cylinder of a 
transient can be raised conventionally, and the Air Fuel Ratio Control precision of a transient can be raised. 
[0028] By the way, at claim 2 and claim 13 which were mentioned above, although the variation of the 
restoration air content in a cylinder to intake valve close timing was predicted based on the prediction throttle 
opening, since operation parameters, such as engine rotational speed, change after predicting the variation of the 
restoration air content in a cylinder based on a prediction throttle opening before [ intake valve close timing ], 
the predictability of the restoration air content in a cylinder falls to a transient under the effect. 
[0029] As this cure, like claim 16, while presuming the current restoration air content in a cylinder based on a 
current throttle opening Predict a future throttle opening and the future restoration air content in a cylinder is 
predicted based on a future throttle opening. The deflection (equivalent to the prediction variation of the 
restoration air content in a cylinder) of said future restoration air content in a cylinder, and said current 
restoration air content in a cylinder It adds to the restoration air content in a base cylinder calculated based on 
the current operation parameter, the final restoration air content in a prediction cylinder is calculated, and you 
may make it calculate fuel oil consumption based on this final restoration air content in a prediction cylinder. If 
it does in this way, the prediction variation of the restoration air content in a cylinder can be calculated with a 
sufficient precision from deflection with the restoration air content in a cylinder by which the present was 
presumed to be a future restoration air content in a cylinder, and the predictability of the restoration air content 
in a cylinder can be improved. Because, the current restoration air content in a presumed cylinder is because 
information, such as the newest engine speed, is taken into consideration. (In claim 2 grade, it is the variation of 
the forecast of the restoration air content in a cylinder, and the value which becomes origin is calculated in the 
past.) 
[0030] 

[Embodiment of the Invention] «operation gestalt (1)» The operation gestalt (1) of this invention is hereafter 
explained based on drawing 1 thru/or drawing 19 . 

[0031] First, based on drawing 1 , the outline configuration of the whole control system of the engine 1 1 which 
is an internal combustion engine is explained. The upstream of the inlet pipe 12 of an engine 1 1 is equipped 
with an air cleaner 13, and the air flow meter 14 which measures an inhalation air content is installed in the 
downstream. The supply current to a heat ray is controlled so that this air flow meter 14 keeps constant the 
temperature gradient of the temperature of a heat ray and the intake-air temperature which the heat ray (not 
shown) and intake-air-temperature sensing element (not shown) which are arranged in the flow of inhalation air 
are built in, and are cooled with inhalation air. According to the heat release of a heat ray which changes 
according to an intake air flow, the supply current to a heat ray changes by this, and the voltage signal according 
to this supply current is outputted as an intake-air-flow signal. 

[0032] A throttle valve 15 is formed in the downstream of this air flow meter 14, and the motors 17 (throttle 
actuator), such as a DC motor, are connected with rotation shaft 15a of this throttle valve 15. And the opening 
(throttle opening) of a throttle valve 15 is controlled by driving force of this motor 17, and this throttle opening 
is detected by the throttle opening sensor 18. 

[0033] In this case, a throttle opening is controlled by driving force of a motor 17, the inhalation air content at 
the time of idle operation is controlled by it also at the time of idle operation, and it carries out feedback control 
so that an engine speed may be made in agreement with target idle rotational speed. The throttle control at the 
time of this idle operation is idle rotational-speed control (ISC). In addition, idle rotational-speed control 
prepares an idle rotational-speed control bulb (ISC bulb) in the bypass path which bypasses a throttle valve 15, 
and it controls the opening of this idle rotational-speed control bulb, and you may make it control a bypass air 
content (inhalation air content at the time of idle operation) at the time of idle operation. 
[0034] On the other hand, the intake-pressure sensor 1 6 which detects an intake pressure is installed in the 
downstream of a throttle valve 15. Moreover, a fuel injection valve 20 is attached in the inlet manifold 19 which 
introduces into each cylinder of an engine 1 1 the inhalation air which passed the throttle valve 15, and the 
ignition plug 21 is attached in the cylinder head of each cylinder of an engine 1 1 . The periphery of the signal 
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rotor 23 attached in the crankshaft 22 of an engine 1 1 is countered, the crank angle sensor 24 is installed, the 
pulse of the engine revolution speed signal Ne outputted from this crank angle sensor 24 is incorporated by the 
electronic control unit (ECU) 25, and an engine speed is detected by the generating frequency of this engine 
revolution speed signal Ne. 

[0035] On the other hand, the amount of treading in of an accelerator pedal 26 (accelerator control input) is 
detected by the accelerator sensor 27, and the voltage signal according to this accelerator control input is 
incorporated by the electronic control unit 25 through A/D converter 28. Moreover, the output of an air flow 
meter 14, the intake-pressure sensor 16, and the various sensors of throttle opening sensor 18 grade is also 
incorporated by the electronic control unit 25 through A/D converter 28. 

[0036] It is performing the various programs for throttle control which this electronic control unit's 25 is 
constituted considering the microcomputer equipped with CPU29, ROM30, and RAM31 grade as a subject, and 
are memorized by ROM30 by CPU29. At the time of the usual throttle control, feedback control of the motor 17 
is carried out by PID control etc. through the motorised circuit 32 according to the opening command value 
(target throttle opening) set up based on the accelerator control input etc. A throttle opening is controlled by 
driving force of this motor 1 7 to an opening command value. In addition, all over the energization way from the 
motorised circuit 32 to a motor 17, the safety circuit 46 which consists of a relay etc. is formed, at the time of 
the abnormalities of an electronic throttle system, this safety circuit 46 operates and the energization to a motor 
1 7 is intercepted. 

[0037] Furthermore, this electronic control unit 25 is performing each routine of drawing 1 0 memorized by 
ROM30 thru/or drawing 1 8 by CPU29. While performing throttle delay control mentioned later, the throttle 
opening of intake valve close timing (definite timing of the restoration air content in a cylinder) is predicted. 
The restoration air content in a cylinder is predicted based on the prediction throttle opening, fiael oil 
consumption is calculated based on the restoration air content in a prediction cylinder, the injection pulse of the 
pulse width according to the result of an operation is outputted to the actuation circuit 45, and the injection time 
(fuel oil consumption) of a fuel injection valve 20 is controlled. 

[0038] The operation approach of the fuel oil consumption by this electronic control unit 25 is explained based 
on drawing 2 thru/or drawing 9 . Drawing 2 is the block diagram showing the outline of the prediction approach 
of throttle delay control and the restoration air content in a cylinder. During engine operation, an accelerator 
control input is detected by the accelerator sensor 27, and an opening command value (target throttle opening) is 
set up by the map etc. according to an accelerator control input etc. with an opening command value operation 
means. This opening command value is predetermined time Tdly with a delay means. It is delayed and is 
outputted to the motorised circuit 32 of an electronic throttle system, time delay Tdly of this opening command 
value it is shown in drawing 3 ~ as — time amount Tinj from the operation timing (prediction timing of the 
restoration air content in a cylinder) of the fuel oil consumption TAU of a certain cylinder to the intake valve 
close timing of the cylinder concerned from — it is set as the time amount (Tdly =Tinj-Tth) which deducted the 
dead time Tth of an electronic throttle system. 

[0039] In this case, the dead time Tth of an electronic throttle system is the time amount Tinj from the operation 
timing (prediction timing of the restoration air content in a cylinder) of fuel oil consumption TAU to intake 
valve close timing, although it does not change even if a throttle actuation rate changes. It becomes so short that 
an engine speed becomes a high revolution. For this reason, at the time of a high revolution, it is the time 
amount Tinj from the operation timing of fuel oil consumption to intake valve close timing. There is a 
possibility of becoming shorter than the dead time Tth of an electronic throttle system. 

[0040] This point is taken into consideration and it is the time amount Tinj from the operation timing of fuel oil 
consumption TAU to intake valve close timing at this operation gestalt (1). When becoming shorter than the 
dead time Tth of an electronic throttle system, it outputs without delaying an opening command value. 
[0041] Opening command value phitotal before making it delayed with a delay means on the other hand It is 
inputted into an electronic throttle model. This electronic throttle model consists of the electronic throttle 
dynamic-model section and variation operation part, as shown in drawing 4 . while this electronic throttle 
dynamic-model section simulates an electronic throttle system response delay property with a second-order lag 
element [omega2/(s2+2 zetaomegas+omega 2)] — the marginal property of the actuation rate of a throttle valve 
1 5 — a rate limiter — simulating — opening command value phitotal before delay from — prediction throttle 
opening thetaf It calculates. Two integral elements (l-/s) of a second-order lag element are rectangle integrals. 
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In addition, in order to simplify data processing, it may replace with a second-order lag element, and a first 
order lag element may be used. 

[0042] Moreover, the variation operation part of an electronic throttle model consists of a derivative element 
(d/dt) and an integral element (integral), is a derivative element (d/dt), and is the sampling time Ts of the output 
(prediction throttle opening) of the electronic throttle dynamic-model section. It asks for the difference of a 
between, and it is integrating an integral element (integral) with this difference, and prediction variation 
deltatheta of a throttle opening is calculated. Under the present circumstances, the time amount which integrates 
an integral element (integral) with difference is the time amount Tinj from the operation timing (prediction 
timing of the restoration air content in a cylinder) of fuel oil consumption TAU to intake valve close timing. It 
is the larger one of the dead time Tth of an electronic throttle system. Thereby, prediction variation deltatheta of 
the throttle opening outputted from variation operation part becomes the prediction variation of the throttle 
opening to intake valve close timing (or after dead-time Tth progress). 

[0043] An electronic throttle model adds prediction variation deltatheta of the throttle opening outputted from 
variation operation part to the current throttle opening theta (output of the throttle opening sensor 1 8), and is 
prediction throttle opening thetaf of intake valve close timing (or after dead-time Tth progress). It asks and is 
this prediction throttle opening thetaf. It outputs to an inhalation-of-air system model. 

[0044] As this inhalation-of-air system model is shown in drawing 5 , it consists of prediction throttle passage 
air content operation part, prediction intake-pressure operation part, and the restoration air content operation 
part in a prediction cylinder, and prediction throttle passage air content operation part considers that throttle 
opening which inhalation air passes is an orifice, and calculates the prediction throttle passage air content Gin 
from a prediction throttle opening etc. moreover, prediction intake-pressure operation part ~ the prediction 
throttle passage air content Gin to prediction intake pressure Pm calculating — the restoration air content 
operation part in a prediction cylinder — prediction intake pressure Pm from — the restoration air content Gcf in 
a prediction cylinder is calculated. Prediction throttle passage air content operation part is expressed with the 
formula of the following orifice. 
[0045] 
[Equation 3] 
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Gin = m • A • • f(Pm/Pa) (1) 

/R-T 

Gin : *nyMoiia^SlS[kg/sec] 

A : ^cyh^BBP^3»BfiB«|[m»] 
Pa : *Mm[Pk) 
Pm : »$tffi[Pa] 

R : ^#5£* 
T : KftiBJSCK] 
A = jtr*(l-cos 2 0) (2) 

f(PWPa)=y_l fr (_2 r )^r ...-(3) 
ic : itmt 

[0046] Here, f (Pm/Pa) is Pm/Pa, in order to simplify data processing, although you may calculate by the top 
type. It is good to compute from the table made into a parameter. The table of f (Pm/Pa) is expressed with a 
curve as shown in drawing 6 when maximum is normalized as 1 . Since this f (Pm/Pa) does not become a 
negative value physically so that clearly from the above-mentioned (3) formula and (4) types, it has been set as f 
(Pm/Pa) = 0 in the example of drawing 6 at the time of Pm/Pa>l . 

[0047] However, when it is made f(Pm/Pa) = 0 at the time of Pm/Pa >1 , there is an inclination which the 
operation value (the throttle passage air content Gin, the prediction intake pressure Pm, and restoration air 
content Gcf in a prediction cylinder) of an inhalation-of-air system model vibrates, and hunting generates at the 
time of heavy load operation to which Pm/Pa is changed in the one neighborhood as shown in drawing 8 . This 
cause is Pm/Pa. Since f (Pm/Pa) is guarded by 0 whenever Pm/Pa becomes one or more on an operation at the 
time of heavy load operation for the rate of change of f (Pm/Pa) to become large in the field of the one 
neighborhood, it is because change of f at the time of heavy load operation (Pm/Pa) becomes irregular. 
[0048] As this cure, with this operation gestalt (1), the table of f (Pm/Pa) is set up, as shown in drawing 7 . 
That is, it is set as an f(Pm/Pa) = negative value at the time of f(Pm/Pa) =0Pm/Pa >1 at the time of f(Pm/Pa) = 
forward value Pm/Pa =1 at the time of Pm/Pa <1 . Thereby, the table of f (Pm/Pa) serves as a symmetrical 
change property which ** reverses bordering on Pm/Pa =1 . 

[0049] If the table of f (Pm/Pa) of a change property like drawing 7 is used, change of f (Pm/Pa) will become 
regular at the time of heavy load operation to which Pm/Pa is changed in the one neighborhood. For this reason, 
by equalizing the operation value (the throttle passage air content Gin, the prediction intake pressure Pm, or 
restoration air content Gcf in a prediction cylinder) of an inhalation-of-air system model, as shown in drawing 
9 , the output (restoration air content Gcf in a prediction cylinder) of the inhalation-of-air system model at the 
time of heavy load operation can be stabilized, and hunting can be prevented. 
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[0050] Intake pressure Pm inputted into prediction throttle passage air content operation part Although the last 
prediction intake pressure Pm (i-1) calculated by prediction intake-pressure operation part is used, you may 
make it use the output of the intake-pressure sensor 16. 

[0051] Moreover, with this operation gestalt (1), although the throttle opening effective sectional area A used 
for the operation of the prediction throttle passage air content Gin may compute the throttle opening theta by 
substituting it for the aforementioned (2) formula, in order to simplify data processing, multiplication value mu- 
A of a flow coefficient mu and the throttle opening effective sectional area A is computed from the table which 
makes a prediction throttle opening a parameter. 

[0052] Next, prediction intake pressure Pm The operation approach of the restoration air content Gcf in a 
prediction cylinder is explained. If the law of conservation of mass is applied to the flow of the inhalation air 
which flows the inhalation-of-air path (henceforth a "throttle down-stream inhalation-of-air path") from a 
throttle valve 15 to the inlet port of an engine 11, the relation expressed with the following (5) types will be 
obtained. 

d/dt-Qm =Gin-Gcf .... (5) 

Here, it is Qm. The air content in a throttle down-stream inhalation-of-air path, and d/dt-Qm A prediction 
throttle passage air content and Gcf of the variation of the air content in a throttle down-stream inhalation-of-air 
path and Gin are the restoration air contents in a prediction cylinder. 

[0053] Moreover, if a gaseous equation of state is applied to a throttle down-stream inhalation-of-air path, the 
relation expressed to the following (6) equations will be obtained. 
Gcf=eta -(Ne/2) -Vc - (Qm/VIM) .... (6) 

eta: Volumetric efficiency Ne : Engine speed Vc : Cylinder-capacity VIM: Content volume of a throttle down- 
stream inhalation-of-air path [0054] Here, volumetric efficiency eta is the engine speed Ne which is a parameter 
with an intake air flow and a correlation in order to change with intake air flows. Intake pressure Pm It is based 
and is set up on a map etc. Pm used here It is the last value Pm of a prediction intake pressure (i-1). 
Eta=f(NeandPm) 

[0055] Moreover, model time constant tauIM of an inhalation-of-air system model is expressed with the 
following (7) types. 

tauIM=2andVIM/(Vc and eta-Ne ) .... (7) 

The following (8) types are drawn from the above-mentioned (5) - (7) type. 
d/dt-Qm =Gin-Qm/tauIM .... (8) 

[0056] The above-mentioned (8) formula is discretized as follows, in order to be able to carry out data 
processing of this with an electronic control unit 25, since it is continuous system. 
{Qm(i)-Qm(i-1)}/Ts -Gin(i) -Qm(i-l)/tauIM .... (9) 
Here, it is Ts. It is the sampling time. 

[0057] If this (9) type is arranged, it is the air content Qm in a throttle down-stream inhalation-of-air path. 

Operation expression is drawn as follows. 

Qm(i) ={Gin (i)-Qm(i-l)/tauIM} and Ts+Qm (i-1) [kg] 

....(10) 

[0058] moreover — if a gaseous equation of state is applied to a throttle down-stream inhalation-of-air path — air 
content Qm in a throttle down-stream inhalation-of-air path from — prediction intake pressure Pm The equation 
to calculate is drawn as follows. 
Pm =Qm and R-T/VIM [Pa] .... (11) 

R: Gas constant T : the above-mentioned (10) formula and (11) types are used for the prediction intake-pressure 
operation part of an intake-air-temperature inhalation-of-air system model, and it is the prediction intake 
pressure Pm. It calculates. 

[0059] The operation expression of the restoration air content Gcf in a prediction cylinder expressed with the 
following (12) types is drawn from the above-mentioned (11) formula and (6) types. 
Gcf=eta-Vc -Pm/(2 and R-T) [kg/rev] .... (12) 

The restoration air content operation part in a prediction cylinder of an inhalation-of-air system model 
calculates the temporary restoration air content Gcf in a prediction cylinder using the above-mentioned (12) 
formula. 

[0060] As shown in drawing 2 , it is inputted into a derivative element (d/dt), the difference between the 
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sampling times ts is called for, and the difference integrates an integral element (integral) with the output 
(temporary restoration air content Gcf in a prediction cylinder) of an inhalation-of-air system model, the reset 
time — time amount Tinj from the operation timing (prediction timing of the restoration air content in a 
cylinder) of fuel oil consumption TAU to intake valve close timing it is . The value with which the integral 
element (integral) was integrated is prediction variation deltaGc of the restoration air content in a cylinder to 
intake valve close timing. It becomes a value considerable the bottom and is this prediction variation deltaGc. It 
adds to the restoration air content Gbase in a base cylinder calculated with the base inhalation-of-air system 
model, and the final restoration air content Gc (restoration air content in a cylinder decided to intake valve close 
timing) in a prediction cylinder is calculated. 

[0061] Next, the operation approach of the restoration air content in a base cylinder is explained. This 
restoration air content in a base cylinder is a current restoration air content in a cylinder calculated based on the 
output (intake air flow) of an air flow meter 14. Therefore, the variation of the restoration air content in a 
cylinder by change of the throttle opening from current to intake valve close timing (definite timing of the 
restoration air content in a cylinder) is not contained in this restoration air content in a base cylinder. Although 
there is an advantage that the operation precision of the restoration air content in a cylinder at the time of a 
stationary is good since the method of calculating the restoration air content in a cylinder generally serves as a 
restoration air content in an intake-air- flow = cylinder from the output of an air flow meter 14 at the time of a 
stationary Since the response delay (for example, response delay according [ the case of the air flow meter 14 of 
a heat type ] to the own heat mass of the sensor section of an air flow meter 14) of an air flow meter 14 exists, a 
transient has the fault that the responsibility of a transient is bad. 

[0062] So, with this operation gestalt (1), in order to raise the responsibility of a transient, the response delay of 
the output of an air flow meter 14 is compensated with a response delay compensator (phase-lead-compensation 
element), the output of this response delay compensator is inputted into a base inhalation-of-air system model, 
and the restoration air content Gbase in a base cylinder which is the output of this base inhalation-of-air system 
model is calculated. The transfer function of this base inhalation-of-air system model is expressed with the 
following first-order-lag type. 

Gbase=l/(l+tauIM-s), output tauIM of the restoration air-content Gdlay response delay compensator in a 
GdlayGbase:base cylinder: Time constant [0063] Time constant tauIM of this base inhalation-of-air system 
model is expressed with a degree type. 
tauIM=2andVIM/(Vc and eta-Ne ) 

content volume Vc of the inhalation-of-air path of the VIMithrottle downstream : cylinder capacity eta: — 
volumetric efficiency Ne : an engine speed ~ engine speed Ne which is a parameter with an intake air flow and 
a correlation here since volumetric efficiency eta changes with intake air flows Based on an intake pressure P 
(output of the intake-pressure sensor 16), it is set up on a map etc. 

[0064] Prediction variation deltaGc of the restoration air content Gbase in a base cylinder calculated with such a 
base inhalation-of-air system model, and the restoration air content in a cylinder calculated from the prediction 
throttle opening etc. It integrates, the final restoration air content Gc (restoration air content in a cylinder 
decided to intake valve close timing) in a prediction cylinder is calculated, and it is this restoration air content 
Gc in a prediction cylinder. Fuel oil consumption is set up according to an engine speed etc. 
[0065] The function of each block of drawing 2 explained above is realized by each routine of drawing 10 
thru/or drawing 18 . Hereafter, the content of processing of each routine is explained to a detail. 
[0066] The main routine of [main routine] drawing 10 is performed a predetermined period after ON of an 
ignition switch. Opening command value phitotal set up according to the accelerator control input etc. when this 
routine was started, the throttle delay control routine of drawing 1 1 mentioned later was performed at step 1 00 
and the execution condition of throttle delay control was satisfied first Predetermined time Tdly Throttle delay 
control to delay is performed. Then, the restoration air content operation routine in a prediction cylinder of 
drawing 12 which progresses to step 200 and is mentioned later is performed, and the restoration air content Gc 
(restoration air content in a cylinder decided to intake valve close timing) in a prediction cylinder is calculated. 
[0067] Then, it progresses to step 300, a basic injection-quantity operation routine (not shown) is performed, 
and it is the restoration air content Gc in a prediction cylinder. Engine speed Ne It responds and is the basic 
injection quantity Tp by a map etc. It calculates. Then, it is the basic injection quantity Tp about various kinds 
of correction factors Kc, such as the fuel correction factor [ as opposed to / perform the injection-quantity 
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amendment routine of drawing 18 which progresses to step 400 and is mentioned later, and / a load effect ] 
Kload (acceleration-and-deceleration correction factor), an air-fiiel ratio feedback correction factor, and a water 
temperature correction factor. Multiplication is carried out and final fuel oil consumption is calculated. 
[0068] The throttle delay control routine of [throttle delay control-routine] drawing 1 1 is a subroutine 
performed at step 100 of the main routine of drawing 10 . If this routine is started, first, it responds to an 
accelerator control input (output of the accelerator sensor 27) etc. at step 101, and is opening command value 
phi total. It sets up. under the present circumstances, opening command value phi total Demand opening phipedal 
according to an accelerator control input Demand opening phiisc by idle rotational-speed control (ISC) etc. — it 
integrates and asks for various demand openings. 

phitotal =phipedal+phiisc [0069] Processing of this step 101 plays a role of an opening command value 
operation means as used in the field of a claim. 

[0070] Then, it progresses to step 102 and judges whether throttle delay control prohibition conditions are 
satisfied, it is that the time of that it is predetermined within a time [ immediately after the time of ** start up or 
start up ] as throttle delay control prohibition conditions here, for example and ** idle operation or an 
accelerator control input is small, that ** automatic transmission is in a neutral condition, etc., if there are 
conditions which correspond at least one of these conditions, throttle delay control prohibition conditions will 
be satisfied, and when other, throttle delay control prohibition conditions become abortive. 
[0071] The time of start up, and immediately after start up, since the engine revolution is unstable, when throttle 
delay control is performed essentially, there is a possibility that engine revolution fluctuation may become still 
larger. Moreover, at the time of idle operation, in order for idle rotational-speed control (ISC) to operate and to 
carry out feedback control of the idle rotational speed, when throttle delay control is performed, there is a 
possibility that idle rotational-speed control may interfere with throttle delay control, and an idle revolution may 
become instability. Moreover, since an operator may perform racing (an engine carrying out a fore blow) when 
an automatic transmission is in a neutral condition, when it is in a neutral condition and throttle delay control is 
performed, the start of an engine speed is overdue at the time of racing activation, and there is a possibility that 
accelerator responsibility and acceleration nature may impress an operator with it being bad. 
[0072] So, with this operation gestalt (1), the adverse effect by throttle delay control is lost by forbidding 
throttle delay control, when it is the operational status in which the adverse effect by throttle delay control 
appears (at the time [ At the at the time of idle operation immediately after start up ] of a neutral). [ At the time 
of start up ] 

[0073] If judged with throttle delay control prohibition conditions being satisfied at step 102, throttle delay 
control will be forbidden, and it will progress to step 103, and will output to the motorised circuit 32, without 
delaying current (newest) opening command value phitotal (i). 

[0074] On the other hand, at step 102, if judged with throttle delay control prohibition conditions being 
abortive, throttle delay control will be carried out as follows by processing after step 104. First, it is opening 
command value phitotal at step 104. Time delay Tdly It determines, under the present circumstances, time delay 
Tdly it is shown in drawing 3 — as — time amount Tinj from the operation timing (prediction timing of the 
restoration air content in a cylinder) of fuel oil consumption TAU to intake valve close timing from — it is set as 
the time amount (Tdly =Tinj-Tth) which deducted the dead time Tth of an electronic throttle system. However, 
time amount Tinj from the operation timing of fuel oil consumption TAU to intake valve close timing When 
becoming shorter than a dead time Tth (in the case of Tinj-Tth<0), it is a time delay Tdly. It is referred to as 0. 
[0075] Then, it progresses to step 105 and is a time delay Tdly. The inner number Cdly of samplings It 
calculates by the degree type. 

Cdly =Tdly / Ts — here — Ts It is the sampling time. 

[0076] Then, it progresses to step 106 and is a time delay Tdly from current. The inner number Cdly of 
samplings Opening command value phitotal (i-Cdly) calculated before is outputted to the motorised circuit 32. 
thereby — opening command value phitotal output timing — time delay Tdly only — you make it delayed 
Processing of these steps 102-106 plays a role of a delay means as used in the field of a claim. 
[0077] The restoration air content operation routine in a prediction cylinder of [restoration air content operation 
routine in prediction cylinder] drawing 12 R> 2 is a subroutine performed at step 200 of the main routine of 
drawing 1 1 , and plays a role of a restoration air content operation means in a prediction cylinder as used in the 
field of a claim. 
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[0078] If this routine is started, first, at step 201, the prediction intake-pressure operation routine of drawing 13 
mentioned later will be performed, and the prediction intake pressure Pm (intake pressure of intake valve close 
timing) will be calculated. Then, it progresses to step 202 and is the prediction intake pressure Pm. It uses and is 
the restoration air content Gcf in a prediction cylinder (i) by the degree type. It calculates. 
Gcf (i) = eta- Vc -Pm/(2 and R-T) [kg/rev] 

eta: Volumetric efficiency Vc : Cylinder-capacity R:gas constant T: Intake-air temperature [0079] Then, it 
progresses to step 203 and judges whether throttle delay control prohibition conditions are satisfied by the same 
approach as step 102 of drawing 1 1 . It progresses to step 204 from step 203, without carrying out throttle delay 
control, when throttle delay control prohibition conditions are satisfied, and is prediction variation deltaGc of 
the restoration air content in a cylinder from the operation timing of fuel oil consumption to intake valve close 
timing. It is set as 0. 

[0080] On the other hand, when throttle delay control prohibition conditions are abortive, throttle delay control 
is carried out by the throttle delay control routine of drawing 1 1 , and it progresses to step 205 from step 203, 
and is the prediction time amount Tinj. The inner number Cp of samplings It calculates by the degree type. 
Cp =Tinj / Ts — here ~ prediction time amount Tinj the time amount from the operation timing of fuel oil 
consumption to intake valve close timing — it is — Ts It is the sampling time. 

[0081] Then, it progresses to step 206 and is prediction variation deltaGc of the restoration air content Gcf in a 
cylinder from the operation timing of fuel oil consumption to intake valve close timing. It calculates by the 
degree type. 

deltaGc = — Gcf (i) -Gcf (i-Cp) — here — Gcf (i) this restoration air content in a prediction cylinder (that is, 
restoration air content in a prediction cylinder of the operation timing of intake valve close timing) — it is — Gcf 
(i-Cp) — current ~ prediction time amount Tinj The inner number Cp of samplings It is the restoration air 
content in a cylinder (that is, restoration air content in a cylinder of the operation timing of fuel oil 
consumption) calculated before. 

[0082] Prediction variation deltaGc It progresses to step 207 after an operation, the restoration air content 
operation routine in a base cylinder (not shown) is performed, and the restoration air content Gbase in a base 
cylinder is calculated. Under the present circumstances, the response delay of the output of an air flow meter 14 
is compensated with a response delay compensator (phase-lead-compensation element), and the restoration air 
content Gbase in a base cylinder is calculated with the following transfer function using the output Gdlay of this 
response delay compensator. 

Gbase=l/and (1+tauIM-s), and Gdlay — here, tauIM is a time constant. Although the operation expression of the 
restoration air content in a base cylinder is expressed with the continuous system in order that a top type may 
simplify explanation, an electronic control unit 25 calculates the restoration air content Gbase in a base cylinder 
using the formula discretized in the top type. 

[0083] Then, prediction variation deltaGc for which progressed to step 208 and this restoration air content 
Gbase in a base cylinder was asked at the above-mentioned step 206 It adds and is the final restoration air 
content Gc in a prediction cylinder. It asks. 

Gc = Gbase+delta Gc [0084] The prediction intake-pressure operation routine of [prediction intake-pressure 
operation routine] drawing 13 is a subroutine performed at step 201 of the restoration air content operation 
routine in a prediction cylinder of drawing 12 . If this routine is started, first, at step 211, the prediction throttle 
passage air content operation routine of drawing 14 mentioned later will be performed, and the prediction 
throttle passage air content Gin will be calculated. Then, the inhalation-of-air system model time constant 
operation routine of drawing 16 R> 6 which progresses to step 212 and is mentioned later is performed, and 
model time constant tauIM of an inhalation-of-air system model is calculated. Then, it progresses to step 213 
and is the air content Qm in a throttle down-stream inhalation-of-air path. It calculates by the degree type. 
Qm(i) ={Gin (i)-Qm(i-l)/tauIM) and Ts+Qm (i-1) — here — Qm (i) — the air content in this throttle down- 
stream inhalation-of-air path, the air content in the throttle down-stream inhalation-of-air path of last time 
[ Qm / (i-1) ], and Ts It is the sampling time. 

[0085] then, the step 214 — progressing — air content Qm in a throttle down-stream inhalation-of-air path from - 
- prediction intake pressure Pm It calculates by the degree type. 

Pm =Qm and R-T/VIM — here, R is [ an intake-air temperature and VIM of a gas constant and T ] the content 
volume of a throttle down-stream inhalation-of-air path. 
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[0086] Then, it is the prediction intake pressure Pm by progressing to step 215 and computing the average of 

this prediction intake pressure Pm (i) and the last prediction intake pressure Pm (i-1). It equalizes. 

Pm (i) = {Pm(i)+Pm (i-l)}/2 [0087] The prediction throttle passage air content operation routine of [prediction 

throttle passage air content operation routine] drawing 14 is a subroutine performed at step 2 1 1 of the prediction 

intake-pressure operation routine of drawing 13 . If this routine is started, the prediction throttle opening 

operation routine of drawing 15 mentioned later is first performed at step 221 , and it is prediction throttle 

opening thetaf of intake valve close timing. It calculates. Then, it progresses to step 222 and atmospheric 

pressure Pa, an intake-air temperature T, and the last prediction intake pressure Pm (i-1) are read. 

[0088] Then, it progresses to step 223 and the prediction throttle passage air content Gin is calculated by the 

degree type. 

[0089] 

[Equation 4] 

Gin = fi • A • ■ Pa • f (Pm/Pa) 



Gin : ^n^hAoffiig^^SCkg/sec] 

Pa : *5aflE[Pa] 

Pm : K^tflECPa] 
R : 

f(Pm/Pa) : &fmPmb*ttPa<Dt£^&*Z®mm 



[0090] Under the present circumstances, mu-A is prediction throttle opening thetaf. Computing from the table 
made into a parameter, f (Pm/Pa) is Pm/Pa. It computes from the table of drawing 7 made into a parameter. 
Intake pressure Pm The last prediction intake pressure Pm (i-1) is used, and, as for atmospheric pressure Pa and 
an intake-air temperature T, the detection value of a sensor is used, respectively. In addition, atmospheric 
pressure Pa A standard atmosphere (fixed value) may be used. 

[0091] The prediction throttle opening operation routine of [prediction throttle opening operation routine] 
drawing 15 R> 5 is a subroutine performed at step 221 of the prediction throttle passage air content operation 
routine of drawing 14 , and plays a role of a prediction throttle opening operation as used in the field of a claim. 
[0092] If this routine is started, first, it responds to an accelerator control input etc. at step 231, and is opening 
command value phitotal. It sets up. under the present circumstances, opening command value phitotal Demand 
opening phipedal according to an accelerator control input Demand opening phiisc by idle rotational-speed 
control (ISC) etc. — it integrates and asks for various demand openings. 

phitotal =phipedal+phiisc [0093] Then, it is opening command value phitotal before delay by the electronic 
throttle dynamic-model section and variation operation part of the electronic throttle model which progresses to 
step 232, progresses to step 233 after reading the present throttle opening theta detected by the throttle opening 
sensor 18, and is shown in drawing 4 . It uses and prediction variation deltatheta of a throttle opening is 
calculated. This prediction variation deltatheta is the time amount Tinj from the operation timing (prediction 
timing of the restoration air content in a cylinder) of fuel oil consumption TAU to intake valve close timing. It 
is the prediction variation of a throttle opening. However, time amount Tinj to intake valve close timing When 
shorter than the dead time Tth of an electronic throttle system, prediction variation deltatheta of the throttle 
opening within this dead time Tth is called for. 

[0094] Then, it progresses to step 234, prediction variation deltatheta is added to the current throttle opening 
theta, and it is prediction throttle opening thetaf. It asks. 

thetaf = theta+delta theta — this prediction throttle opening thetaf It is the prediction throttle opening of intake 
http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_webcgi_ejje 1 0/1 7/2005 



JP,2002-201998,A [DETAILED DESCRIPTION] 



Page 14 of 17 



valve close timing (or after dead-time Tth progress). 

[0095] The inhalation-of-air system model time constant operation routine of [inhalation-of-air system model 
time constant operation routine] drawing 16 R> 6 is a subroutine performed at step 212 of the prediction intake- 
pressure operation routine of drawing 13 . If this routine is started, first, at step 241, the volumetric-efficiency 
operation routine of drawing 17 mentioned later will be performed, and volumetric efficiency eta will be 
calculated. Then, it progresses to step 242 and model time constant tauIM is calculated by the degree type. 
tauIM=2andVIM/(eta-Ne [ VC and ]/60) 

Here, VIM is the content volume (fixed value) of a throttle down-stream inhalation-of-air path, and VC. 
Cylinder capacity (fixed value) and Ne It is an engine speed (rpm). 

[0096] The volumetric-efficiency operation routine of [volumetric-efficiency operation routine] drawing 1 7 is a 
subroutine performed at step 241 of the inhalation-of-air system model time constant operation routine of 
drawing 16 . Starting of this routine first reads the last intake pressure Pm (i-1) and the last atmospheric 
pressure Pa, an intake-air temperature T, an engine speed Ne, valve timing WT, and the cooling water 
temperature THW at step 151 . Then, basic volumetric-efficiency etar progress to step 152, search the 
volumetric-efficiency map which makes a parameter Pm/Pa, an engine speed Ne, and valve timing WT, and 
corresponding to the current engine operation condition It calculates and is this basic volumetric-efficiency etar. 
It amends with the correction value according to the cooling water temperature THW, and volumetric efficiency 
eta is searched for. 

[0097] The injection-quantity amendment routine of [injection-quantity amendment routine] drawing 1 8 is a 
subroutine performed at step 400 of the main routine of drawing 10 , and plays a role of a fuel-oil-consumption 
operation means as used in the field of [ routine / (not shown) / basic injection-quantity operation ] a claim. 
[0098] a ****** [ whether it is a load effect (fluctuation of the restoration air content in a cylinder) by 
accelerator actuation at step 401 , and that for example, an accelerator control input is beyond the set point first 
when this routine is started ] — or it judges by whether the variation of an accelerator control input is beyond the 
set point. It is the small value Kl about the fuel correction factor [ as opposed to / if judged with the load effect 
by accelerator actuation / progress to step 402 and / a load effect (fluctuation of the restoration air content in a 
cylinder) ] Kload. It sets up. By the operation approach of the restoration air content in a cylinder of this 
operation gestalt (1), since this reason can predict the load effect (fluctuation of the restoration air content in a 
cylinder) by accelerator actuation with a sufficient precision, it is because amendment to fuel oil consumption 
can be lessened. 

[0099] It is the large value K2 about the fuel correction factor [ as opposed to / when judged with it not being a 
load effect by accelerator actuation (in for example, the case of the load effect by the time of shifting an 
automatic transmission from a neutral range to a drive range, power steering, a brake, an air-conditioner, etc.) / 
progress to step 403 and / a load effect on the other hand ] Kload. It sets up. Since this reason cannot predict the 
load effect by factors other than accelerator actuation from an accelerator control input, it is because it is 
desirable to make [ many ] amendment to fuel oil consumption to the load effect by factors other than 
accelerator actuation. 

[0100] After determining the fuel correction factor [ as opposed to a load effect as mentioned above ] Kload by 
step 402 or 403, it progresses to step 404 and various kinds of fuel correction factors Kc to factors other than a 
load effect (for example, an air- fuel ratio feedback correction factor, a water temperature correction factor, a 
study correction factor, etc.) are calculated, and they are the basic injection quantity Tp and the fuel correction 
factors Kload and Kc at the following step 405. And invalid injection time Tv It uses and the final fuel oil 
consumption (injection pulse width) TAU is calculated by the degree type. 

TAU=Tp xKloadxKc+Tv[0101] An example of the behavior of the prediction throttle opening calculated by 
each routine explained above and the restoration air content in a prediction cylinder is shown in the timing 
diagram of drawing 19 . under engine operation — an accelerator control input etc. ~ responding — opening 
command value phitotal setting up — this opening command value phitotal output timing ~ time delay Tdly 
only — you make it delayed under the present circumstances, time delay Tdly it is shown in drawing 3 — as — 
time amount Tinj from the operation timing (prediction timing of the restoration air content in a cylinder) of 
fuel oil consumption TAU to intake valve close timing from — it is set as the time amount (Tdly =Tinj-Tth) 
which deducted the dead time Tth of an electronic throttle system. However, time amount Tinj from the 
operation timing of fuel oil consumption TAU to intake valve close timing When becoming shorter than the 
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dead time Tth of an electronic throttle system (in the case of Tinj-Tth<0), it is a time delay Tdly. It is referred to 
asO. 

[0102] Opening command value phi total before delay It is based, prediction variation deltatheta of a throttle 
opening is calculated with the electronic throttle model of drawing 4 R> 4, this prediction variation deltatheta is 
added to the present throttle opening theta (output of the throttle opening sensor 1 8), and it is prediction throttle 
opening thetaf of intake valve close timing (or after dead-time Tth progress). It asks. And this prediction throttle 
opening thetaf It uses, the temporary restoration air content Gcf in a prediction cylinder is calculated with the 
inhalation-of-air system model of drawing 5 , differential / integral processing of this is carried out, and it is 
prediction variation deltaGc of the restoration air content in a cylinder to intake valve close timing. It calculates. 
This prediction variation deltaGc It adds to the restoration air content Gbase in a base cylinder calculated with 
the base inhalation-of-air system model, and the final restoration air content Gc (restoration air content in a 
cylinder decided to intake valve close timing) in a prediction cylinder is calculated. Thereby, it is the restoration 
air content Gc in a cylinder. It can become possible to predict with a sufficient precision, and the Air Fuel Ratio 
Control precision of a transient can be raised. 

[0103] «operation gestalt (2)» With the above-mentioned operation gestalt (1), although the delay means 
(refer to drawing 2 ) was made to perform throttle delay control, in the operation gestalt (2) of this invention 
shown in drawing 20 , a delay means is excluded, and throttle delay control is not performed, but a throttle 
opening is predicted using the dead time Tth of an electronic throttle system. 

[0104] With this operation gestalt (2), it outputs to the motorised circuit 32, without delaying the opening 
command value set up based on the accelerator control input etc. with the opening command value operation 
means. By the same approach as said operation gestalt (1), and with an electronic throttle model Based on an 
opening command value and the present throttle opening (output of the throttle opening sensor 1 8), the throttle 
opening of intake valve close timing (or after dead- time Tth progress) is predicted. With an inhalation-of-air 
system model (configuration of drawing 5 ) The temporary restoration air content in a prediction cylinder is 
calculated from this prediction throttle opening, differential / integral processing of this is carried out, and the 
prediction variation of the restoration air content in a cylinder to intake valve close timing (or after dead-time 
Tth progress) is calculated. And it adds to the restoration air content in a base cylinder which calculated this 
prediction variation with the base inhalation-of-air system model, and the final restoration air content in a 
prediction cylinder is calculated. 

[0105] Also in this operation gestalt (2) explained above, it can become possible to predict a throttle opening 
using the dead time Tth of an electronic throttle system, and to predict the restoration air content in a cylinder 
with a sufficient precision from the prediction throttle opening, and the Air Fuel Ratio Control precision of a 
transient can be raised. 

[0106] «operation gestalt (3)» Although the above-mentioned operation gestalt (1) and (2) apply this 
invention to an engine with an electronic throttle system, the operation gestalt (3) shown in drawing 21 applies 
this invention to the engine with a mechanical-cable-type throttle system which interlocks a throttle opening 
with accelerator actuation mechanically. 

[0107] With this operation gestalt (3), since an accelerator control input and a real throttle opening interlock 
mechanically and there is no response delay of a throttle, the opening command value operation means, the 
delay means, and the electronic throttle model are not prepared. Although the prediction throttle opening was 
inputted into the inhalation-of-air system model, with this operation gestalt (3), a current throttle opening 
(output of the throttle opening sensor 1 8) is inputted into an inhalation-of-air system model said operation 
gestalt (1) and (2). The configuration of this inhalation-of-air system model is substantially [ as said operation 
gestalt (1) ] the same, calculates the temporary restoration air content in a prediction cylinder from a current 
throttle opening, carries out differential / integral processing of this, and calculates the prediction variation of 
the restoration air content in a cylinder to intake valve close timing (or after predetermined period progress). 
And it adds to the restoration air content in a base cylinder which calculated this prediction variation with the 
base inhalation-of-air system model, and the final restoration air content in a prediction cylinder is calculated. 
[0108] Also in the case of a mechanical-cable-type throttle system, if it does in this way, the operation precision 
of the restoration air content in a cylinder can be raised conventionally, and the Air Fuel Ratio Control precision 
of a transient can be raised. 

[0109] «operation gestalt (4)» in said operation gestalt (1) and (2) As shown in drawing 2 and drawing 20 , 
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based on an opening command value and the present throttle opening (output of the throttle opening sensor 18), 
the throttle opening of intake valve close timing (or after dead-time Tth progress) is predicted with an electronic 
throttle model. With an inhalation-of-air system model Calculate the temporary restoration air content in a 
prediction cylinder from this prediction throttle opening, and differential / integral processing of this is carried 
out. After calculating the prediction variation of the restoration air content in a cylinder to intake valve close 
timing (or after dead-time Tth progress), he adds to the restoration air content in a base cylinder which 
calculated this prediction variation with the base inhalation-of-air system model, and is trying to calculate the 
final restoration air content in a prediction cylinder. 

[0110] on the other hand, with the operation gestalt (4) of this invention shown in drawing 22 thru/or drawing 
25 With an electronic throttle model, the throttle opening of intake valve close timing (or after dead-time Tth 
progress) is predicted based on an opening command value and a current throttle opening (output of the throttle 
opening sensor 18). While calculating the future restoration air content in a cylinder (temporary restoration air 
content in a prediction cylinder) from this prediction throttle opening, with an inhalation-of-air system model 
with an inhalation-of-air system model Based on a current throttle opening (output of the throttle opening 
sensor 1 8), the current restoration air content in a cylinder is calculated. The deflection (equivalent to the 
prediction variation of the restoration air content in a cylinder) of the future restoration air content in a cylinder, 
and the current restoration air content in a cylinder It adds to the restoration air content in a base cylinder 
calculated with the base inhalation-of-air system model, and he calculates the final restoration air content in a 
prediction cylinder, and is trying to calculate fuel oil consumption based on this final restoration air content in a 
prediction cylinder. 

[0111] Except for the point of performing the same main routine as drawing 10 , and performing the restoration 
air content operation routine in a prediction cylinder of drawing 23 at step 200 also by this operation gestalt (4), 
it is the same as said operation gestalt (1). 

[0112] By the restoration air content operation routine in a prediction cylinder of drawing 23 , the present 
restoration air content presumption routine in a cylinder of drawing 24 mentioned later is performed by step 
500, and it is based on the present throttle opening theta (output of the throttle opening sensor 1 8) with an 
inhalation-of-air system model, and is the present restoration air content Gest in a cylinder. It calculates. 
[01 13] then, the future restoration air content operation routine in a cylinder of drawing 25 which progresses to 
step 600 and is mentioned later — performing — an electronic throttle model — an opening command value and 
the present throttle opening theta — being based — throttle opening thetaf of intake valve close timing (or after 
dead-time Tth progress) predicting — an inhalation-of-air system model ~ this prediction throttle opening thetaf 
from — the future restoration air content Gcf (temporary restoration air content in a prediction cylinder) in a 
cylinder is calculated. 

[0114] Then, restoration air content Gest in a cylinder progress, the future restoration air content Gcf in a 
cylinder, and current to step 700 after progressing to step 600 and calculating the restoration air content Gbase 
in a base cylinder by the same approach as said operation gestalt (1) Deflection (equivalent to the prediction 
variation of the restoration air content in a cylinder) is added to the restoration air content Gbase in a base 
cylinder, and it is the final restoration air content Gc in a prediction cylinder. It asks. 
Gc =Gbase+(Gcf-Gest ) 

[01 15] In the present restoration air content presumption routine in a cylinder of drawing 24 , the present 
throttle opening theta is first read at step 501 , and they are atmospheric pressure Pa, an intake-air temperature T, 
and an intake pressure Pm at the following step 502. It reads. Under the present circumstances, intake pressure 
Pm The detection value of the intake-pressure sensor 16 may be used, or the last value of the prediction intake 
pressure calculated at step 601 of drawing 25 mentioned later may be used. 

[01 16] Then, it progresses to step 503 and the present throttle passage air content Gin is calculated by the same 
approach as the routine of drawing 14 explained with said operation gestalt (1). Then, it progresses to step 504 
and model time constant tauIM of an inhalation-of-air system model is calculated by the same approach as the 
routine of drawing 16 explained with said operation gestalt (1). 

[01 17] Then, it is the air content Qm in a throttle down-stream inhalation-of-air path with the same approach as 
step 213 of drawing 13 which progressed to step 505 and was explained with said operation gestalt (1). It 
calculates by the degree type. 

Qm(i) ={Gin (i)-Qm(i-l)/tauIM} and Ts+Qm (i-1) — here — Qm (i) — the air content in this throttle down- 
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stream inhalation-of-air path, the air content in the throttle down-stream inhalation-of-air path of last time 
[ Qm / (i-1) ], and Ts It is the sampling time. 

[0118] then, the step 506 — progressing — air content Qm in a throttle down-stream inhalation-of-air path from - 
- current intake pressure Pm It calculates by the degree type. 

Pm =Qm and R-T/VIM — here, R is [ an intake-air temperature and VIM of a gas constant and T ] the content 
volume of a throttle down-stream inhalation-of-air path. 

[01 19] Then, it is an intake pressure Pm by progressing to step 507 and computing the average of this intake 
pressure Pm (i) and the last intake pressure Pm (i-1). It equalizes. 

Pm (i) = {Pm(i)+Pm (i-l)}/2 [0120] Then, it progresses to step 508 and is this intake pressure Pm. It uses and is 
the current restoration air content Gest in a cylinder by the degree type. It calculates. 
Gest = eta-Vc -Pm/(2 and R-T) 

eta: Volumetric efficiency Vc : Cylinder-capacity Rigas constant T: Intake-air temperature [0121] On the other 
hand, in the future restoration air content operation routine in a cylinder of drawing 25 , the prediction intake 
pressure Pm (intake pressure of intake valve close timing) is first calculated at step 601 by the same processing 
as the prediction intake-pressure operation routine of drawing 13 explained with said operation gestalt (1). 
Then, it progresses to step 602 and is the prediction intake pressure Pm. It uses and the future restoration air 
content Gcf (restoration air content in a cylinder of intake valve close timing) in a cylinder is calculated by the 
degree type. 

Gcf=eta-Vc -Pm/(2 and R-T) 

[0122] With this operation gestalt (4) explained above, while presuming the current restoration air content in a 
cylinder based on a current throttle opening Since a future throttle opening is predicted, the future restoration air 
content in a cylinder is predicted and the prediction variation of the restoration air content in a cylinder was 
calculated from the deflection of said future restoration air content in a cylinder, and said current restoration air 
content in a cylinder From said operation gestalt (1) and (2), the prediction variation of the restoration air 
content in a cylinder can be calculated with a sufficient precision, and the predictability of the restoration air 
content in a cylinder can be improved. 

[0123] «operation gestalt (5)» Although the inhalation-of-air system model which calculates the restoration 
air content in a cylinder from a throttle opening was used in above-mentioned operation gestalt (1) - (4), with 
the operation gestalt (5) of this invention, air content change sets the time constant of this inhalation-of-air 
system model as a small value also for a actual twist using the inhalation-of-air system model which calculates 
the restoration air content in a cylinder from the output (intake air flow) of an air flow meter 14 (intake-air- flow 
detection means) so that it may appear at an early stage. 

[0124] With this operation gestalt (5), the time constant of an inhalation-of-air system model is set as a small 
value instead of establishing a means to predict the restoration air content in a cylinder. If it does in this way, 
the effectiveness as predicting the future restoration air content in a cylinder that a actual twist also has it since 
change of the restoration air content in a cylinder calculated with the inhalation-of-air system model appears at 
an early stage will be acquired. [ same ] Thereby, the operation precision of the restoration air content in a 
cylinder of a transient can be raised conventionally, and the Air Fuel Ratio Control precision of a transient can 
be raised. 

[Translation done.] 
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* NOTICES * 
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damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 
2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 
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Pm =Qm • R * T/V XM [Pa] (1 l) 

R : «flc^» 30^ [0 0 5 9 ] ±15 (11) iCt (6) 



T : ttftflUE 

Kas^^/Ko^sjK^EEffltKSfPtt, ±15 do) at 

(1 1) *«rfflv^T, ^'M^EEPra-£rS®[-r5o ♦ 
G c f • J£j^jl£. C2-V.R 
!K«K^^O^fflt»rt3c*^ftft»»ffl5«i. ±15 (1 
2) Sfct/Hl^T, (El5a^«^«^JW.cf 
o e 

[0 0 6 0] i2 (c^-TJ: p ic, Kags^T^u-ofcHA 
«Eo^9JfBrtjl5«£«4Gcf) f*. fifc&KSS (d/d 
t) KAASiv f-^y >^b#K1 t s BJ<aSi$M**«> 40 
&*t, (J) -C«#S*l3. ^rOfJS 

^*-COB#MTinj -CfcS. ffiftSX (J) T'M^L/i 

AGc Sr-^-^.SSS^7 f /HCj:t)}S{JUfc-<— J*fUrt 
XttSftSGbaselCftlJK LT, »l*Wfc^a!l«rt56tf[£ 

^tt) so 



(12) x-cssns^nnrtXttSMGcfotRiFsc: 

-T-)~ [kg/re v] (1 2) 

[0 0 6 1 ] W:, ^—^iart3E«3E«ft<oiR»*feS: 

*14^fcti;^ (RAffiJMK*) (ca-5t^T«*L^affi 

co^^att (w*.tf»i5s;<oj.T7ci>— ^ i Acorns 

[0 0 6 2] -tc-e, *H«a^fii ( i ) "Cli. 
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Gbase = 1 / ( 1 + t rM • s ) • Gdlay 
Gbase : Xffiftftt&ZE^vk 
Gdlay : ttVgftttffKXattM 
r im 

[0 0 6 3] :<o^- j*K«^7vu£>b#£$: t XM (4, 

r i M =2 • Vjm/ (Vc • t) • Ne ) 
Vim : *° y h/uTSElWoa«iBB<OrtSat 

vc >» >ymn 

Ne : ^^^VSlciUg 
[0 0 6 4] ;*!SMX^^A4Cj:9ft* 

[0 0 6 5] KJiKMUfcH2 0*^oy^©fl|l|gtt. 
HI 0 75MH11 8 0#yu— ^>{CctorHm^ix5o 

[0066] c^^r ^/u— m i o<7>^ -< wu— 
x4„ -r^^^^a >-^^ *y^(D*>-mz.mj£fflmT~n 

oor% ^yi-r^m 1 1 <t>* □ ^ H^aswa/u-f-y 

total SrBf^B*BBTdly iSJI S * 5 * ° h/HSSWJ» 

[0067] ^7^30 o\zm^. m*m&i 

.lei**, aSS-J-SHl 8 0Pftfitfl:MijE^-^>Srll?T 
U ftffiKttK:#r«tt»4t£«ftKload GtoKiSffijE 

«0#fill^4ilijE«Sf[Kc4-S*PaWffiTp «C»»LT* BO 
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[0 0 6 8] [^ny b/^M&MW^-^ V] 01 
^fy/l 0 OT-Hfr S*t5*^— ^V"C*>3. 

~C\ RCtt«|^ total $:S:^-r5 0 M«Jf4*» 
total (4, T^-feyuttf^ftcffilSfcK^BBS^peda 
1 t. T-r K*«!fifcXMff (ISC) »dJ:5Sf*H« 
10 <t>isc «?<D**ftSf*M««:8[KLU-C*«)S. 
<f> total = <f> pedal + <J> isc 

[0 0 6 9] :^7y/i o i <Dj&w.frV?Ptt&tf:<offi. 

[0 0 7 0] 77^1 O 2(Cti^ 7Dyh 

r-r K/^5iteB#xj4r^-fe'/^f^fi^/hSv^i, <3)i 
20 boftft© 5 %o^f i o-etRa t5ftw 5 fcji 

[0071] ttn&*»ttlbiE«r4» ^Wtci^v^^iHi 

r-r K/hi(eb#j4, r>< hvnaeitflMSP (i s 
-rztztb. v^&mmm&fto k/h§he 

— h7/u^cfigOct^(4, jSte#*5u— (xy^y 

[0 0 7 2] -r^X% *HJS^ffi (1) T-14. 

[0 0 7 3] t> L, ^77/10 2^ ^nyh/ugjg 

(**r) o|»a»«tt[* total (i)tiBif **-ric^~# 
®EE/jIeIK3 2 JctB^i"-5« 
[0 0 7 4] — ^7^/1 0 2T\ 5/ h/UgJi 

WSrStJSi"*. *-f, ^7^/l04t', MflEIS-frffi* 
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total (DSiP#fWTdly COR. iliPSIH 

Tdly li, i3iOTfJ:9i:, J»f*®WfiTAU<Df»W 

ft'VU/H*^ ^>^*T*OB*IBTinj 
h/u^fAW *MUVIfl T t h £ S L 3 1 V * fcB^ (Tdly 
= Tinj -Tth) {C|3!^$fL5o fiU JB» B Jttt*TA 

B*BTinj ^ftHkB#B0TthJ:9 tH< (Tinj 
-Tth<0O4§^) J*, i5fiB*MTdly £0<tr£o 
[0 0 7 5] ^f^y 1 0 StCiH^ S3£P$M 10 

Tdly rtOlfi'/y^ftCdly «fc 0 Mfci" 

So 

Cdly =Tdly /Ts 

£r.-c\ Ts 

[0 0 7 6] 77 7/ 1 0 6lC*i*« %tte£V 

fciB3&*lfflTdly rtOI^/y r^ftCdly ffllwtRJt L 
total (i-Cdly) *5Mft[E]K3 2 

mz>-rz>o zti{cx<9^ mmm^m* total otitj^^-r 

/1 0 2-1 0 6 04t&3S|AS«»m* #<E>ffiE8t?^ 9 7-< U 20 

[0 0 7 7] [^««rt3E*ffi**W»y^— HI 1 
y^-fyo^f ^^2 0 0T*Hft ^ixS-y-Z/i—^:/ 

[0 0 7 8] ^-VsftSjEMSixSi:, £^777/ 

2 o i t% «a&-t-aia 1 3 (D^mmn.E.mn^—^^^ 

Mil U ^»J»«J£ P m (««'</u/B3 

»E) £rfi3?i-5o dtf>&, 777/2021:1*11^ ^ 30 

Gcf(i) $r^i"5 0 

Gcf (i) = 77 • Vc • Pm / ( 2 • R • T) [k r/ 
rev] 

vc : ->y ^y^ia 

R : Salfc^gfc 
T : 5&«iaffi 

[0 0 7 9] 777/2 0 3fwii^ x il 10 

^•^102^1:^ 70 y hJ^mmMmm± 40 

^fW^iLT^5^^l»JSt6 0 707 h/uJSffi 
Igj&Srfr-f l-*7 7/2 0 3^^fy7'2 0 4 Jdifi^ 

[0 0 8 0] 7 a y h/UiB]&Mff$tJkAf^3dS7dc 

ft0>»<frl±, HI l l^ny h/i'8BE»]UM — f-^lcj; 

?>77 7/2 0 SlCii^ ^ffl]D*raTinj rtOlflx/y 50 
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Cp =Tinj /Ts 

[0 0 8 1 ] wO^, 77^/2 0 6UliS^ ifW* 

ISIrtasaSEftftGcfo^as^kftAGc 

AGc =Gcf(i) -Gcf(i-Cp) 

::t\ Gcf(i) tt, (o£ 

3zm&%&) -CfctK Gcf(i-Cp)tt, affi^Ot^SJH* 

[0 0 8 2] WJffttAGc <D$tg*£. 77 7/2 0 
-T) *3*fTU 7«rt^«ffi**Gbase£StJ*i- 

O^3tLffi«S^<0l±S^GdlayS:^V>T, ft^fiSH 
SKJ: 0^— 7«rt^*£«*Gbase&jRJ|^rS. 
Gbase = 1/ (l + ti M - s) • Gdlay 

fc*tfflv^-C-<"-7lBrt3c9[SS«a:GbaseS:at»i-S. 

[0 0 8 3] "tCD'te^ 777/2 0 8 fwig-^, 
^«rt**ffi*aGbaselC. JiiE77 7/2 0 6t'^ 

Gc =Gbase+ A Gc 
[0 0 8 4] [fflI®IE^I^-f>] HI 1 3 (O^m 

^/u— fy^777/2 o l T'Htr^ns-y-//^— ^> 
^>a*fiB&SJx5i. tf^T7/2 l 

1 T% ^ili-r^HI 1 4 <7>^a'J7 a y 
/u-fy$rlffl, fffl^ny K/Wiii^^^aGin^r?S 
W--fZo ^(om. 777/2 l 2 tdifg^ ^xBi-^HJi 

6 cOK^^^7^B$^ScrSS/^"-^>^:^tTL.. 
^7/1^^7/^$^$: r iM^r^Si'-So ^^^s 77 7 
/2 1 3l^ii^, 7n-/ h/UTg[55^ii^rt^^^SQ 

Qm(i)= {Gin(i) -Qm(i-l)/r i«) • Ts +Qm(i- 
1) 

ww-C, Qm(i)(l^lHl07n y h^TS£!8«i§ifeF*i<0ffi 

^s. ts ji-^^/y >tn#r$X'3bZ> 0 

[0 0 8 5] wCO?g, 77 7/2 1 4 7D;yh 
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Pm =Qra - R • T/V IM 

-~t\ Ritm.frm&. Tttaafttaa. VxMft^n^h 

[0 0 8 6] 77^2 1 5(Cii^ % ^(O^ 

WK*EEPm(i)fc«riao^«»*JEPm(i-l)i0^tt 

Pm (i) = {Pm(i) + Pm(i-1)} /2 

[0087] a * vj^^%.m^n^-^- 10 

Pa 



Gin = 
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Si 3 ©TMRJIUEWIM'-- fy^fyy2 1 1 

*a y H/uBaSSIIM*— *-v«:3«tU K«^</^G3^ 

tf(mE!«>^«»«EP»(i-l)fcK*atr. 
[0 0 8 8] CLCO^, s/^2 2 3lca*. f»^n 

[0 0 8 9] 
f(Pm/Pa) 



Gin 

A 
Pa 
Pm 



/R-T 

*j<UE[Pa] 
«^£E[Pa] 



f(Pm/Pa) : RSUEPxn^*«JBEPaOtt;-e**5«>affl 



[0090] cois, m • a^. fn^oy h/Mwao 

f fc^9>— * ^Ud»feHttU f (Pm/ 

Pa ) li, Pm /Pa > — * £ i"-5El 7 CO^T— 7 

/u^e>»ffl^5. »«E-Pin fi, l»E^M*»EPin 30 
(i-l)riS/Bv>t>*x, *ft/EPatK«tifiST«i, *ti*fti 

[0 0 9 U fmb* o y b;i>mmmn/\ — ^v] la 1 

5 o^W* a ^ h/I'MftStll/l'— 7- Vft, il4 0^$] 
j*d^ h/ui§ig^^fi:^g[/u~^>co^^^y 2 2 1T' 

[0 0 9 2] */i — Six* 

23 1T% r^-t^|(kff*««C«i:i:THia«*«*tota 40 
1 ^l^^-T5o Z<0^ DRttfl|«fl[« total It, T^ir 
/u»^iB:»wjcSDfeS*BBffi0 pedal fc % T-f K^EHES 
ft»J» (ISC) KJ:5R3fcHIS*isc eo«*ftK* 

<i> total =0 pedal + <f> isc 

[0 0 9 3] :Og x 3 2JCit^, ^nyh 

tt*&A>*:». 3 3 Kil^. EUlo^m^ 

^bftWSffltlc.toraSHlF^BBlSffi^fi* total 50 



^4MMb*A 0 (i. »PWJttTAU<0«*^^ ^ >^ 

Tinj li^Dy K/U->^^A<OjSlflfcB*WTthJ:9 
[0 0 9 4] wO^, -7^2 3 4 {ZM^. &1£<D* 

0f = 0 + A 0 

^ (XttJKft«pB3TthSia^) (Dfid^oyh/Hatt' 

[0 0 9 5] ^ V] [U 1 

6 <0««S*tvub5£R«JM'-- Bi 3 co^fflij 

y^2 4 i-c\ tg^-rsdi 7 (omMZhmmw^- 

r i„= 2 • V,„/ (V c • tj - Ne /6 0) 



(14) 

25 

M) . Vc teisyyymm mizm . Ne ^^^^^ 

(r pm) T*£>5 0 

[0096] ifrmzhmmn^-?' ^] in ^^fg^ 

"C\ ftu0OB5^JEPm(i-l). ^EPa . SS^U^T, 

Tk^THW^rH^iAtfo wCOm. ^X^^15 2tCig 10 
Pro /Pa , ^V^EfwSgNe % x</U^^>f ^ 

lifeffljEfl|[-CttjE-bT»«a*i) Sr*«>-5. 

[0097] [p*w*ffiE/u— ^ ei 1 8 (Dm&$m$ 

-f > (B*iW) fc*^»W»*Ol8IBTU^5«»P|l 
Wi«**fttLto8llfc*fct. 20 
[0 0 9 8] ^■^SjgjftSix^i:, ttT^y-f 

4 0 IT. T**/U»f«C («Srt*«ffl« 

o 2 Mit^K A'ffisn («rt36«ffi«*oa») 
sat** *mmxm (i) ©«rt3E*2E*a[*^«»*tt 

[0 0 9 9] — T?-tM&mc «fcS:&^£b-c&^ 

»*fijE«ftKloadSr**V^K2 CKSt5c reDS& 
/^akfP**»?>tt^»J-C**v^«>. T**/u«m^G> 40 

[0 10 0] W-bOj: 5 (CUT, ^fi'^O 2X^:4 
0 3T\ *WSBB^»i-5«»«E«ftKload*»36L 

it'&. *"ry?4 0 4{ciS^. £^Sba*<^HU:*f 
i-5*a«)«3|sfllB]E««Kc («*.!* SS««:7 -f - K'< 
y^»3E«», *iMB3Eff», *W»jEflWM» 
U ^^f^7 , 4 0 5T\ I*«WtTp x ^trfilE 
^SfcKload. Kc Rt/JSaiBWBSBflTv £ffl<^T, ffirl& 



tfrBH 2002-201998 
26 

TAU = Tp XKloadxKc + Tv 
[0 10 1] »±Ift?*Lfc*A^^>i£J:oT*l(Lfc 

Srl2l 1 9<D?s( Aft- MC^LT^5 e ^Vv^S'Sc 
*tt. r^-fe/u»f^**»ca:DTH«*S*«[* total £ 
BU6U wOM«»««* total Oti^*^ S V^Sril 
JPSMTdly *fttSS£-fr£. iiJP^MTdly 

^^A^MI*:B#r^Tth$:^U5!^^cB^ (Tdly =Tin 
j -Tth) KKJESn*. ffl.U H8R«WiTAUO» 

inj S^ny h/U£x;*^AOlSttB#PflTthJ: 9 fc 
fi<tt5»fr (Tinj -Tth<0O»-fr) li. SKftM 
Tdly £0 tt^o 
[0 10 2] SSIttonatt^I^ total (CS^^TIH 

lAGc ^r^c— ^Si^^^^UCct^^^Ufc^— 
tfiS^fi) $r3ft465o CtxirJ: D , tfil^?c^^^:SGc 

[0103] «n*£^sg (2) » ±ia^ffi^fii ( 1 ) T' 

fWSrfT^ J:5i:Lfe^ HI 2 0 \z^1-^$£W(DM1&&B 
(2) T*(i, ^ WfSS:«v^ h/HgJESi] 

[oio4] *mmmm (2) r-(i. iw«ts*«[«w^ 
lt, iRrexttfl^ffi (1) iB*o*s-e, «^oy 

«^<yu^H^-f (Xtt«tt«FlBTthttiBlfc) 
a 5/ h/uDBflEfe^WL* KSI^f^ (EI5 0«fiJt) 
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[0 10 5] «±RW Lfc#j|flBBai (2) IC*5V*T 

e>«rt**3**4:ttft a < TBtSC £ a* pTIE t ft 

[0106] «H«S^9§ (3) » ±EXfi»B ( 1 ) . 

(2) tt. *«IBS:fI-T-^o y h;Uv/^f^Wt^x> 

(3) tt. hA^BI!flE«:T^-fe/u»f1?^««W(Cig 
[0 10 7] *5li£jgffi (3) -Ctt % T^^J^m^&t 

iS^ffi (i) . (2) -m, ffists^e^/nc^aj^B y 
h^Ha*A*u^*«. *nm&m (3) -en. sa^ 

fffiSSHlfiflgffi (1) irUKWlcra 3i.?£<0* d 

[0 10 81 rcoj: a«^Dj/ hA->* 

[0109] atibiMi (4) » mexttBS (1) . 
( 2 ) -e«t, m 2 % m 2 o \ztMt x o i-. ° y h 

✓^^ph*^ 5^ (xiiaitkB^BaTthaiaa) <^>^° 40 

*uy h/uM^fefio^»JfBrt3EJKS«**i!t»L.. 

(xtt««kB*WTthSia«) *-cotairt5C«ffi«fto^ 

[0110] ^ni^HT, 12 2 2 SlC^-f* 

&w<o%vmWi (4) T'rt, tf^py h/^r/wa so 
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fi) fcfltJH-S-*, SJM^f/HaoT, 

hAHH* h/HMK-t^lM 80tB^) ICS 

[0 111] *JtlK)MI (4) Tt, ilO^D^ 
y/i/- 7->«:JtfTU * J ?y'?2 0 0-CB2 3<O^WfB 
rt3E*ffiK**»^-*-^«r*tT-f-aA«rl»v>-C. sola 
JUKJMI (1) £P]CT*£>5« 

[01 12] in 2 z^mm^m^n.m^n^—^^ 

[oi 13] ^^^^6 o o\zm*h. 

(Xtt«ftB#MTthaiSS) 0* a y h/H58ffi 0f 

[0 1 14] ^co^. 0 OlciSI^., mjfSHJS 

JF^flg (1) trai^^ifeT^— ^!Brt*«S«ftGbase 

2s«*Gcfts. j aE<oi8rt**ffi**Gest t<omm m 

«*GbasetC*D»LT. *»Wft^fflffirt3E«ffi«fi:G 
c ^^^^ e 

Gc =Gbase+ (Gcf— Gest ) 

[oi is] i2 4<D^tE<om^m^n.&i^^—^ 

«ta2T&lMS«EEPm SrR^iatr. KEP 
iE-r5IH2 5^^7^y^6 0 1 '?ffig--tZ>=f'®m%iE<n 
[oi 16] rois, ^T'^S 0 3(CiI^- % ffJE^SS 

mm (i) -ciftMLfcia 1 4^^— ^^tmmnjjm 

X\ T&tECOTsU y h^iffi)S^^SGin$:®Si"5o 

ft, ^77/5 0 41:1^ HMfS^ffiJKfg (1) T'lfcW 
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[0 117] ^GD&, ^x^^5 0 5(Cil^. fiufS|y£ 
Km (1) T^U:il307T77'2 13^iO 

Qm(i) = (Gin(i) - Qm(i-l)/ r rM } ■ Ts + Qm(i- 
1) 

Qm(i-l)ttWSO^o y h/UT«EKftii!»rtOffi 10 
Sift. Ts tel^^y >'^B»IB'CfcS. 

[0118] ^"ryZfS 0 6(cii^ x y h 

yuT«K»ftiiKrt^ffi«*Qm a»e>W5E0>JMU£Pin £ 
8t*«^J: 9 

Pm =Qm • R • T/V IM 

Cw-C% R Ttt«M&flE* V IM l^n 7 h 

[0 119] rofgL 0 7triH^ v ^Ih]cd.® 

^/E Pra(i) k *ils]<0»SU£ Pm(i-l) <h *> mi" 
Z^kX\ ffiSUEPm Sr¥^bi"5o 20 
Pm (i) = { Pm(i) + Pra(i-1)} /2 

[0120] :og, ^fyy5 0 8Ci^ ^o®m 

t 

Gest = t, • Vc • Pm / (2 • R • T) 

Vc : v-y is?&m 

R : ^ffrT&Sfc 
T : 

[0121] -£\ El 2 5«o#*o«rt3E«ffi»*«l(: 30 

(1) -CRO^LfeEll 3 0^M8S«ffija«/u— ^^il^ 

;/OK«JE) rcoSL y:/6 0 2fcii 

^WgftJKEEPm Srffli^T, ftzttwi: 9 #*ofifi*j;fc; 
*2£«ftGcf (KK^A^H^-f ^ ^fcoWiftftm&fx 

Gcf=7j • Vc • Pm / (2 • R • T) 
[0 12 2] £X_hKMLfc*3«S®l8 (4) T'fi, 

(i) . (2) J:9fc«rt**ffi»*^«ffi{b*t» 
«A<*«>5wt35S-ct, «rt5lE*ffi«fto^B!]»ffi«: 

[0 12 3] «HASJgftg (5) » ±KMSJBffi (1) - 
(4) XI*. y h/u|JBffi^fejKrt3E«[ffi«ffl[SriRa 
i-SK***^tr«ir^fc^ ^P^colliS^Sg (5) 50 
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ft wb* jib £ 9 1 -p-^ t-mn * «t 5 K'h £ v ^m^&m 

[0 12 4] *%ft}gtt (5) Tfi. f§il*33fcg£§lfi;S: 

^n5fc26, wjfcoiBrt^ffiftftSr^JHi-sotiai: 

J»*a s »e>*L*. mtcj: 9 > fi«P*<af«rt*c*£«ft 

[m©cofiH**iftp^] 

nan *&w<Dni&xm (i) 

[12] ^mxm (i) ^MjP3L=yFo«i^ 

[§3] *n y h/uSJSftiJWtfSJWrtiftJKSSCft 
.M*a? M'BME) oam*^ 5 >^*KWi-5*-f * 

[[24] ^^^ay hyU^fVutrTM-.^o y^lftEl 

[ins] p^^^-^/w^r^-r^n y twm 

[126] f (Pro /Pa ) (D^-^U^tK^tC^i-ID 
(^(75 1) 

[12 7] f (Pm /Pa ) (DT—7;U$:m.&&){^-t\g 

[128] WAWaM5B#KIE 6 © f (Pra/Pa ) (^f- 
^/u^ffli/^T^Ufc^S , Jfart56*S^ftGcf^»j5r 

[[2 9] iSft?Sai(EB#lcia 7 <D f (Pm/Pa ) <7}^r — 

-jfr&m^xmn ufc^»jfBrt**2SftftGcfo*»*: 
[i2io] ^xD^m^iti^^-ty^— 

[[2 1 1 ] * n y h/HS5£»JW^— ^><7}&S<Ow?t*x£r 
[12 1 2] ffirt3C«ffi«tft«»^— ^>^ffl.aoacttS: 

imi 3] ^»j««EEffl»/u—^>o«ia^8sns:^-r 

[[2 14] ^jSJKa/EiKJJ/^—^vo&a^SEixSr^i" 
[[2 15] ^?sy^ a y h/uiffiiSggfift»i»^— ^^<o«i 
[[216] ^Rn^^M^WM^— s ?l'<nVm.<OWL 

[[2 18] PAWftfflilE/u— ^>0*Dfa<^SEnS:*-**7 a 
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1 1 • (MUKMH) , 12-a«f, 

IK 1 9—BWS^=*— h\ 2 0-»IM*. 2 5 

7— T^-tr/Hr^iJ-. 3 *S2ffib[H]8go 

[(Hi ] 




Pm/Pa 



[020] mmmm (2) ^a^-^h^K^ 

\mzw mmm o> otwiasy 

[022] mmm (4) o^sijffljj.ny hcoM^^ 

vtWkm io 
[ni2 3] mmm (4) <D^ay«rt£*ffi«*»#/u 

[024] mmm (4) omaEoiart5t«6^«a«i^ 



(18) 
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[12] 







\' 






d 
dt 



mm 



Raws*- 



[0 3] 



gJB^WTdly 
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